Size Effects on Transport Properties in Topological Anderson Insulators 
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We study the size effects on the transport properties in topoiogical Anderson insulators by means 
of the Landauer-Biittiker formahsm combined with the nonequihbrium Green function method. 
Conductances calculated for serval different widths of the nanoribbons reveal that there is no longer 
quantized plateaus for narrow nanoribbons. The local spin-resolved current distribution demon- 
strates that the edge states on the two sides can be coupled, leading to enhancement of backscatter- 
ing as the width of the nanoribbon decreases, thus destroying the perfect quantization phenomena 
in the topological Anderson insulator. We also show that the main contribution to the nonquantized 
conductance also comes from edge states. Experiment proposals on topological Anderson insulator 
are discussed finally. 
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Recently, the study of the topological insulator (TI) 
has triggered great research interestsi"— . It was first 
proposed both from topological band theory and topo- 
logical field theory, and soon experimentally realized in 
HgTe/CdTe quantum wells (QWs)^!^. Although it is in- 
sulating in the bulk, conduction is allowed on the bound- 
ary due to the presence of edge states. These edge states 
are protected by the time reversal symmetry. Each edge 
state is accompanied by its time reversal partner, and 
the number of pairs is an odd number for topological 
nontrivial phase, which leads to an odd-integer quantized 
conductance on each edge. However these stories are re- 
stricted to where the TI is semi-infinite. When two edges 
are getting close to each other, finite size effect^i plays 
an important role. The overlap between edge states from 
opposite edges opens an energy gap so that properties of 
TI can be subtly modified. 

The concept of TI can be generalized into many other 
insulators, such as topological Mott insulator—, topolog- 
ical superconductor—^—. Its generalization in Anderson 
insulator (AI) is also addressed in disordered systenti^"— 
recently. It is believed that the disorder-induced edge 
states are topologically protected and play a central role. 
In this sense such insulator is named as topological AI 
(TAI). Thus, a question arises naturally that whether 
this phenomenon will also be emergent in the narrow 
nanoribbon with finite size. It is interesting to explore 
the role edge state plays in the finite size system. 

In the present paper, we studied size effects on TAI in 
detail. We choose to study the typical topological band 
insulator formed by HgTe/CdTc QW*i, where spin-orbit 
coupling is encoded. TAI phase is addressed at certain 
random strength when the Fermi surface is at the bulk 
conduction band, where a quantized conductance is ob- 
served. We show that the conductance is no longer quan- 
tized plateaus for narrow nanoribbons. To understand 
clearly the physics of this picture, we presented the lo- 
cal spin-resolved current distribution in the disordered 
bar, which demonstrated that the main contribution to 



the nonquantized conductance comes from edge states. 
However, due to the truncation of the coherence length 
between two edges by the finite sample size, disorder can 
induce the interedge scattering. As a result, the TAI 
phase will be suppressed. Through detailed size depen- 
dence study, we found that by decreasing the width of 
the nanoribbon, the coupling between edge states will 
lead to exponential enhancement of the backscattering 
which destroys the TAI phase eventually. 

As a starting point, we consider a HgTe/CdTe QWs 
narrow nanoribbon. The low energy electron states 
are approximately described by an effective four band 
Hamiltoniani: 
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where /i(k) = e(k) + d(k) • ct, k = {kx,ky) is the two- 
dimensional wave vector, a = (iTa:, o^j,, (t^) are Pauli 
matrices. Up to the lowest order of k, d(k) = 
{Akr,,Aky,M - Bk'^), and e(k) = C - DP, where the 
parameters A, B. C and D depend on the thickness of 
HgTc/CdTc QWs. k) is nothing but the time rever- 
sal counterpart of /i(k) so that time reversal symmetry 
is respected. The Hamiltonian is obtained by reducing 
the eight-band Kane model to the reduced Hilbert space 
|£;i,l/2 >, \Hl,3/2 >, |£;i,-l/2 >, and \Hl,-3/2 >. 
Mass M is an important parameter describing the en- 
ergy gap between conduction and valence band. The 
Hamiltonian with M > describes a conventional band 
insulator, while M < corresponds to the TI. In the 
present study, tight-binding lattice model is used, so that 
the above effective model is compactified by substitu- 
tions ki — > -sm{kia), and kf — > ^(1 — cos(fcia)), where 
i ^ x,y, and a is the lattice constant. The width and 
length of the nanoribbon under study are Ly and re- 
spectively. Meanwhile, we introduced disorders through 
random on-site energy with a uniform distribution within 
[~W/2, W/2], with the disorder strength W. 
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FIG. 1: (Color online) The conductance G (a)-(c) vs disorder 
strength W for different Fermi energy Ej. The widths Ly of 
the nanoribbon are (a) Ly = fOOnm, (b) Ly — 200nm, and 
(c) Ly — 300nm. 
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FIG. 2: (Color online) The conductance G vs bar width Ly at 
disorder strength W = 115meV. Different lengths are chosen 
with I/a; = 1000nm (black solid line), 2500nm (red dash line), 
4000nm (green dot line), and 5000nm (blue short dash line), 
while the symbols ■, •, * and i are corresponding to the real- 
istic results of exact diagonalization, respectively. The inset 
shows the linear fitting of ln(2-G) by picking a typical length 
Lx = SOOOnm. 



To calculate the transport properties, we applied the 
method by nonequilibrium Green function-- . A small ex- 
ternal bias V = Vl — Vr is applied longitudinally between 
the two terminals. The local current between neighbor- 
ing sites i and j is calculated by the formulapi^"— : 

J.^, = n^a.JpGlf^^,JEf)]{VL - Vn)] (2) 

where V^l(h) describes the voltages at the lcad-L(R). 
G'"'{Ef) = G^TlC^ is electron correlation function with 
line width function r^(^) = — E2(_r)]i ^^.d the re- 

tarded Green functions CiEf) = [G''{Ef)]'' = l/[Ef - 
Hcen — 5]£ — S^j , with Hcen the Hamiltonian in the cen- 
tral region. The local spin-resolved current Jj^j between 
neighboring sites i and j with spin index a can also be 
calculated from Eq. ^ without summing over spin index 
a. The net current flowing through the central region 
is calculated by summing index i for local currents Ji^i+x 
over any cross-section. After obtaining the current J^, 
the linear conductance is given by G = hniy^o dJL/dV. 
In addition, the linear conductance can be directly ob- 
tained by G = rrpiGTflG"]. 

In the following numerical calculations, we used the 
realistic material parameters of the HgTe/CdTe QWa^: 
A = 364.5mcV-nm, B = -GSGmeV-nm^, G = OmcV, 
D ^ -512mcV-nm2, and M = -lOmcV. The length of 
the nanoribbon = lOOOnm, and the lattice constant 
a = 5nm. Since this model is only valid for small k, we 
set the Fermi energy small around the F point. In the 
presence of disorder, the conductance is averaged over up 
to 400 random configurations except for Fig. |4]where 800 
random configurations are used for each data. 

We first studied the conductance G versus disorder 
strength W for different Fermi energy Ef and the width 



Ly of the nanoribbon, as shown in FigIT] As long as the 
Fermi surface Ef = O.OmeV (see FigH]) lies inside the 
bulk gap, the conductance maintains a quantized value 
2e^//i for different widths in the clean limit. The conduc- 
tance remains this quantized value in a broad range of the 
disorder. When the disorder further increases, electronic 
states become localized, and the conductance decreases 
to zero rapidly. Such observation agrees well with the 
previous result that the quantum spin Hall eSec^^— 
is robust against weak disorder and independent of the 
width of the nanoribbon. It's well known that a finite 
size of the Hall bar opens a gap on the edge due to over- 
lap between edge states^. Thus the edge states are no 
longer Dirac particles. Our result shows that although 
the edge states gain mass in nanoribbon, the topology 
is still preserved. Backscattering is still forbidden due to 
time reversal symmetry according to a brilliant argument 
in Ref. m. 

When the Fermi surface lies in the valence band of 
the bulk such &s Ef = — 14mcV shown in Fig. [TJ sys- 
tem is metallic in the absence of disorder. Once the dis- 
order is turned on, electrons become localized, and the 
conductance dives to zero directly (see Fig. [T]). This 
result exhibits the conventional Anderson localization 
phenomenon. On the other hand, when the Fermi en- 
ergy Ef is raised up into the bulk conduction band, say 
Ef — 14meV, and the width is sufficiently large [see 
Fig. [TJb) and (c)], conductance G decreases gradually 
when the strength of the disorder increases. Comparing 
to the Fig. [T]Jb) and (c), it is apparent that the "dip" 
in conductance (which occurs at ^ QQmeV) becomes 
less and less pronounced as the nanoribbon width Ly in- 
creases. Because the increasing of the nanoribbon width 
Ly, the number of conducting channels in the bulk is in- 
creased. The more conducting channels there are, the 
larger the conductance will be. In addition, the relation 
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FIG. 4: (Color online) (a) the conductance G vs disorder 
strength W for different lengths of the nanoribbon. (b) the 
localization length ^ vs disorder strength W . The Fermi 
energy Ef = 14.0meV and the width of the nanoribbon is 
Ly — lOOnm. 



FIG. 3: (Color online) The distribution of local spin-resolved 
current with disorder strength (a) W = 75.0meV, (b) W — 
90.0meV , and (c) W = llSmeV. The Fermi energy is Ef — 
14.0meV, and the bar size is Lx ~ 200a, and Ly — 20a. The 
vector length is proportional to the square root of the current 
value. 



between conductance G and the nanoribbon width Ly 
obeys: GL^ ~ (^cLy, where CTc is the conductivity, which 
is width independent. Beyond a certain onset of W, the 
conductance turns back and increases to an approximate 
quantized value (2e^//i). G maintains this value for a 
certain range of W before eventually decreases. Concern- 
ing the spin degeneracy, each spin component contributes 
one conductance quanta only. This odd-integer valued 
conductance plateau indicates that the the system is in 
a topologically distinct phase on this stage. More impor- 
tantly, different from the pure TI, the quantized value is 
induced by disorder. This topological nontrivial phase is 
exactly the TAI phase. However, one should note that 
when the width is as small as Ly = lOOnm, as shown 
in Fig. IHa), although the plateau is still present, it is 
no longer integer quantized, but irrationally valued in 
unit of e^//i. At the same time, the conductance plateau 
evolves into a hump structure. 

To further understand the finite size effect on TAI, the 
relation between hump peak value and bar size is inves- 
tigated in detail, as shown in Fig. [5J With the increase 
of width, the plateau conductance G increases as well, 
and finally saturates to 2e'^/h. After numerous calcula- 
tion, we found that the quantity (2 — G) behaves as an 
exponential function of the width Ly, which is quite un- 



conventional. In the case of normal metal, conductance 
is proportional to the width due to the increase of chan- 
nels. This exponential behavior of conductance suggests 
the conductance hump is contributed also by the edge 
channel of the Hall bar, as the interaction between states 
localized on opposite edges decays exponentially with re- 
spect to the width. Due to previous studies by Groth et 
alr^, finite disorder would reverse the sign of mass term 
of the BHZ model^, leading to a robust edge state. The 
above result adds new insights to this understanding in 
narrow nanoribbon structure. 

In order to get a better insight into the microscopic ori- 
gin of the conductance variations as shown in Fig. [IJa), 
the local spin-resolved current distributions are plotted 
in Fig. [3] Due to the time reversal symmetry, we only 
consider the spin- up component, the influence of spin- 
down component can be directly obtained by time re- 
versal symmetry. Here, the local spin polarized current 

on site i is defined as jj = jj^i+x + -^i-^i+y- When 
W = 75.0meV, the whole sample behaves similarly as 
a conventional metal with large backscattering. The 
vortex-like circulation pattern can be easily seen, which 
reflects that the scattering direction are determined by 
the chirality of spin orbital interaction. As to the trans- 
port property, this is an AI as the conductance decreases 
very rapidly as an exponential form with respect to the 
bar length L^ [see Fig. 4(b)]. As the disorder strength 
increases gradually, current on one edge is greatly sup- 
pressed, while enhanced on the other edge. Meanwhile, 
the bulk of the sample becomes more and more insu- 
lating. All these signals indicates a formation of sta- 
ble edge states. At last, when the disorder is raised 
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up to = 115.0meV and conductance reaches peak 
value, bulk states are extremely disordered and insulat- 
ing. Edge states on the upper edge survives and appears 
to be very robust against disorder. However, due to the 
finitcncss of the nanoribbon width, small backscattering 
is still present by hopping from the upper edge to the 
lower one. The reduction of conductance is apparently 
proportional to this hopping probability which is an ex- 
ponential function of the width. 

In Fig. IH^a), the conductance G as a function of disor- 
der strength W with different length of the nanorib- 
bon are plotted. The dip feature is clearer for a large 
because of the increasing probability of the backscatter- 
ing between the two edges. It is interesting to point out 
that the main shape of the conductance hump at around 
the disorder strength W = 115mcV changes little with 
the bar length, and the system in that region exhibits 
the nonquantized version of TAI behavior. The corre- 
sponding localization length ^ = —lirriL^^oaLx /^t^{G) 
as a function of disorder strength W is shown in Fig. 
|31Ib). We notice that the localization length increases 
dramatically near the conductance hump region, which 
is the finite size version of the picture that reason- 
able disorder scattering can drive an insulating system 
into a TAI phase so that the system shows quantized 
conductancei^ii^. On the other hand, <^ almost vanishes 
near the dip region. From the theoretical viewpoint, the 
sharp difference of longitudinal localization length in dif- 
ferent disorder regions plays the central central role for all 



experimental observable effects, both for bulk and finite 
size version of TAI. 

In summary, the size effects on the transport properties 
in TAI are studied in this paper. We found the conduc- 
tance plateau deviates from its quantized value exponen- 
tially as the nanoribbon width decreases. Such behavior 
is originated from the mutual interaction between edge 
states, as is clearly evidenced by local spin-resolved cur- 
rent distributions in the disordered bar. Furthermore, 
we found the conductance is generally decaying exponen- 
tially as the system length increases. However, in the TAI 
region, the longitudinal localization length is extremely 
long, which is found up to hundreds of microns. This is 
quite a good feature that makes the TAI phase with a 
narrow nanoribbon experimentally distinguishable. 

It is rcportcdSi that topological insulator can be re- 
alized in optical lattices, where the disorder can be in- 
troduced by an optical laser speckle potentia P^i^^ . By 
tuning the laser intensity, it's quite promising to observe 
the TAI and the finite size effect. However in contrast 
to the uniform random disorder studied in this work, the 
disorder in optical laser speckle is correlated. The effect 
of this difference is unknown up to now. Therefore, it 
should be an interesting problem for future studies. 
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